Antiferromagnetic spin-1/2 chains in (NO)Cu(N03)3: a microscopic study 
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We report on the microscopic model of the recently synthesized one-dimensional quantum magnet 
(NO)Cu(N03)3. Applying density functional theory band structure calculations, we obtain a leading 
antiferromagnetic exchange coupling J ~ 200 K, which runs via NO3 groups forming spin chains 
along the b direction. Much weaker couplings J' ~ 2 K link the chains into layers in a non-frustrated 
manner. Our calculations do not support the earlier conjecture on an anisotropic frustrated square 
lattice physics in (NO)Cu(N03)s. In contrast, the model of uniform spin chains leads to a remarkably 
good fit of the experimental magnetic susceptibility data, although the low-temperature features of 
the intrinsic magnetic susceptibility measured by electron spin resonance might call for extension of 
the model. We outline possible experiments to observe the suggested long-range magnetic ordering 
in (NO)Cu(NC>3)3 and briefly compare this compound to other spin-^ uniform-chain systems. 

PACS numbers: 75.30.Et, 71.20.Ps, 75.10.Pq, 75.10.Jm 



Quantum magnets give an exciting opportunity to 
observe unusual ground states and to establish unex- 
pected connections between theoretical models and real 
systems*^ Quantum spin chains are in the focus of 
numerous recent studies and show peculiar excitation 
spectra 3 along with the promising effect of the ballis- 
tic heat transport; 4 - Among the spin-chain models, the 
properties of the uniform spin chain are now well under- 
stood theoretically and extensively verified experimen- 
tally via versatile investigations for a range of model 
compounds £&r— The crucial present-day task is to ex- 
tend these results to other systems with different dimen- 
sionality, different lattice topologies, and, consequently, 
different physics. One of the possible approaches to this 
challenging problem is to explore compounds with un- 
usual chemical features that can lead to peculiar crystal 
structures and spin lattices. Yet, the deduction of the 
correct spin model for a complex crystal structure will 
often require a microscopic study to provide quantitative 
estimates of the individual exchange couplings. 

The (NO)Cu(NC>3)3 compound^ is one of the recent 
examples for a low-dimensional magnet with a special 
chemical feature, the nitrosonium [NO] + cation that 
forms a mixed salt with the magnetic spin-i Cu +2 . The 
peculiar crystal structure (Fig. [IJ is formed by chains 
of isolated CUO4 plaquettes running along the b direc- 
tion. One type of the triangular [NO3] - nitrate anions 
links the plaquettes within a chain, thus connecting to 
two neighboring plaquettes. The nitrate groups of the 
second type are connected to one plaquette only. The 
chains stack along the a and c directions, whereas the 
[NO] + cations are found between the chains. The crystal 
symmetry is monoclinic (space group P2\jm). 

An experimental study^ of (NO)Cu(N03)3 evidenced 
low-dimensional spin correlations and proposed a two- 
dimensional (2D) Nersesyan-Tsvelik models which is 
better known as an anisotropic frustrated square 
latticed Based on phcnomenological arguments, in par- 
ticular, i) the almost temperature-independent values of 
the (/-factor and ii) the lack of sharp anomalies in the 
specific heat, the absence of long-range ordering (LRO) 



at least down to 2 K is proposed^ At the same time, the 
experimental magnetic susceptibility evidences that the 
leading magnetic exchange coupling J (along the struc- 
tural chains, Fig.[I| exceeds 150 K. To reconcile the large 
coupling with the absence of LRO, the authors of Ref. 9 
suggest that the interchain couplings J' and Ji in the 
(NO)Cu(N03)3 structure show an exactly 2:1 ratio (see 
Fig. [T]), thus leading to strong frustration that inhibits 
LRO. 

Since individual exchange couplings are not directly 
measurable, it is generally difficult to judge whether a 
specific system is magnetically frustrated or not. For 
instance, such discussion for the spin-1/2 diamond chain 
system azurite CU3 (003)2(011)2 (Ref. 12) is still not set- 
tled: while inelastic neutron scattering data favor non- 
frustrated magnetism,— band structure calculations sug- 
gest a frustrated model ) 14 i 15 whereas thermodynamical 
data can be satisfactorily described by both modelSi 15 i 16 
Even more illustrative is the recent evidence 17 of a 
non-frustrated spin model in (CuCl)LaNb207, initially 
proposed to imply magnetism of the frustrated square 
latticed 

Moreover, the conjecture on the exact J 2 ■ J' — 1 : 2 
ratio in (NO)Cu(N03)3 is based on two non-trivial as- 
sumptions: i) the J 2 and J' couplings are running exclu- 
sively via NO groups; ii) the energies of these couplings 
are proportional to the number of bridging NO units 
(two for J' and one for J 2)^ Regarding the complex- 
ity of exchange interactions in general, such assumptions 
should be supported by a microscopic verification. Band 
structure calculations are known as a reliable and accu- 
rate tool to investigate magnetism on the microscopic 
levelfi^— In particular, this method has been success- 
fully applied to frustrated square lattice systems^ and 
to a variety of Cu +2 compounds! 7 - ' 20 ^ 4 Therefore, we per- 
form band structure calculations and evaluate individual 
exchange couplings in (NO)Cu(N03)3. We find that the 
simple counting of the bridging [NO] + groups is an inap- 
propriate approach, because it does not regard the differ- 
ent geometry of the superexchange pathways^ Our cal- 
culations describe (NO)Cu(N03)3 as an essentially one- 



dimensional (ID) and non- frustrated system. This con- 
flicting finding calls for reconsideration of the available 
experimental data. 

The scalar-relativistic density functional theory (DFT) 
calculations were performed using the full-potential 
f plo9 . 00-33 code i 26 ^ 27 For the local density approxima- 
tion (LDA), the Perdew-Wang parametrizatior* 2 ^ of the 
exchange-correlation potential was chosen. LDA calcu- 
lations were done on a converged mesh of 1920 fc-points 
(588 points in the irreducible wedge). 

LDA is known to fail describing the insulating prop- 
erties of cuprates. Nevertheless, it provides reliable in- 
formation on the relevant orbitals and leading antifcr- 
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FIG. 1. (Color online) Crystal structure (top and bottom) 
and the spin model (middle) of (NO)Cu(N0 3 ) 3 . The neigh- 
boring CuCU plaquettes are connected via NO3 triangles and 
form chains along b (top). The chains are well separated by 
[NO] + cations (bottom). In the top panel, the nearly over- 
lapping NO3 triangles he in different planes and remain dis- 
connected (see also the bottom panel). 
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FIG. 2. (Color online) Top: the valence band of 

(NO)Cu(NOs)3. The Fermi level ef is at zero energy. Bottom 
left: the band structure of the two-band dpa complex at ef 
and the fit using the Wannier functions technique. Bottom 
right: the orbital-resolved density of states for the antibond- 
ing band. 



romagnetic (AFM) exchange couplings. Among others, 
magnetic excitations have the smallest energy, and essen- 
tials of magnetism are concealed in the close vicinity of 
the Fermi level ep. A sharp peak of NO states appearing 
at 1 eV above ef is a peculiar feature of (NO)Cu(NC>3)3 
related to the antibonding 7r*-states of the [NO] + cation 
(two nearly degenerate orbitals for each NO group). In 
other respects, the valence band of (NO)Cu(NOs)3 com- 
prises features typical for cuprates: it has a width of 
about 5 eV and consists predominantly of Cu 3d and O 
2p states (Fig. [2]). The well-separated density of states 
for the antibonding Cu-0 bands at ef has two distinct 
maxima (van Hove singularities), characteristic of a fD 
behavior. Assuming the simplest nearest-neighbor chain 
scenario, the width W of the antibonding band readily 
yields the leading hopping term t = W/i «180 meV. 

To account for all the possible exchange couplings in 
(NO)Cu(N03)3, we consider the valence bands in more 
detail. The two-band complex at ef is formed by a- 
overlapping Cu 2>d x 2_ y 2 and O 2p x>y orbitals (Fig. [2] 
bottom right). The strong hybridization of these or- 
bitals allows to treat them within an effective one-orbital 
model. The band structure (Fig. (3J bottom left) exhibits 
the predominant dispersion along X-S and Y—T. This 
corresponds to the crystallographic b direction and sup- 
ports the proposed ID scenario. To evaluate the hopping 
terms, we ht the valence bands using Wannier functions 
(WF's)22 based on Cu 3d x 2_ y 2 states (Fig.[3]). 3 £ The per- 
fect fit to the LDA band structure (Fig. [21 bottom) justi- 
fies the WF procedure. This way, we obtain t = 150 meV 
for the leading nearest-neighbor intrachain hopping and 




FIG. 3. (Color online) Fragment of the Heisenberg chain. The 
Wannier function for the Cu 3d x 2_ y 2 orbital is shown. 



a small non- frustrated inter-chain hopping t' ~ 17 meV 
(Fig. [U middle). Other hoppings are below 10 meV. In 
particular, the previously proposed ti (see Fig. [T]) ap- 
peared to be as small as 2 meV, disfavoring the model 
with the frustrated interchain couplings. The couplings 
J and J' form layers, whereas the leading hopping in the 
perpendicular direction is t± = 6 meV. 

In contrast to the apparent insulating behavior of 
(NO)Cu(NOs)3 evidenced by the blue color of the 
crystals^ LDA yields a metallic ground state. This short- 
coming of LDA originates from the well-known underes- 
timate of strong electronic correlations, intrinsic for the 
3c? 9 configuration of the magnetic Cu +2 ions. To restore 
the insulating ground state, we add the missing part 
of correlations in two ways: (i) by mapping the LDA 
band structure onto a Hubbard model; (ii) by treating 
the correlations in a mean-field way via the LSDA+J7 
approach with the around-mean-ficld double-counting- 
correction scheme. 

Adopting the mapping procedure, we transfer the lead- 
ing hoppings onto a Hubbard model with the effective 
on-site Coulomb repulsion U e g. Here, the low-energy ex- 
citations can be described by a Heisenberg model, since 
the t <C U e g condition and the half-filling regime are both 
well justified for undoped cuprates^I Assuming a typi- 
cal value of U c s — 4.5 eV) 24 ' 32 > 33 we readily obtain the 
AFM part of the exchange integrals using the expression 
of second-order perturbation theory J AFM = 4i 2 /C7 e ff. 
This way, we find J AFM = 230 K and J' AFM = 3 K. The 
frustrating coupling J AFM ~ 0.04 K is negligible. The 
interlayer coupling is J± — 0.4 K. 

As an alternative approach, we apply the LSDA+J7 
scheme to evaluate J and J', since long-range terms 
are negligible as demonstrated above. The similarity of 
(NO)Cu(N0 3 ) 3 to other Cu+ 2 oxides allows to adopt 
the typical values of the Coulomb repulsion and ex- 
change parameters Ud — 6.5 ± 1 eV and Jd = 1 eV, 
respectively 34 This parameter set yields accurate esti- 
mates of individual exchange couplings for related Cu +2 
compounds* 24 ' 32 ! 33 For (NO)Cu(N03)3, we obtain the in- 



sulating ground state (band gap E g — 1.7 eV) with ex- 
change couplings J = 200 =F 50 K and J' below 1 K. 
Therefore, the model and the LSDA+C/ approaches con- 
sistently describe (NO)Cu(NOs)3 as a ID system with 
the leading exchange coupling of about 200 K and the 
interchain coupling below 3 K. 

The interchain coupling J' likely runs via the NO 
groups, as evidenced by small tails of it* NO molecular 
orbitals in the Cu-based WFs (Fig. |3J). Nevertheless, the 
hoppings depend on the mutual orientation of the WFs, 
hence a simple counting of the bridging NO units neglects 
a basic ingredient of the superexchange mechanism. In 
contrast, our extensive DFT calculations suggest J 2 <C J' 
and do not support the earlier conjecture on the exact 
J' I J% = 2:1 ratio£ In conflict with Ref. [|| we find 
that (NO)Cu(N03)3 is an essentially non-frustrated ID 
spin system. It is rather similar to other ID Cu +2 com- 
pounds with CUO4 plaquettes separated by non-magnetic 
groups. The J value of 150 — 250 K is typical for the Cu- 
O-O-Cu superexchange, e.g., in the uniform-spin-chain 
compounds M 2 Cu(P0 4 ) 2 (M = Sr, Ba)^i 

In the following, we reconsider the experimental data 
for (NO)Cu(N0 3 ) 3 in light of the non-frustrated ID 
spin model suggested by the DFT calculations. We 
first discuss magnetic susceptibility. The susceptibility 
curves were computed via quantum Monte-Carlo (QMC) 
method using the loop 3 ^ and worm algorithms, imple- 
mented in the ALPS simulation package. 36 We performed 
simulations for finite lattices with periodic boundary con- 
ditions. The typical lattice size was N = 60 — 100 for 
ID models and N = 1500 - 2000 for the model of cou- 
pled spin chains. Calculations for lattices of different size 
showed negligible finite-size effects for the temperature 
range considered. Regarding the experimental data, we 
first discuss the temperature dependence of the electron 
spin resonance (ESR) intensity that can be taken as a 
direct measure of %£ 

The ESR data show a broad maximum at T^ ax ~ 
100 K. In the uniform-spin-chain model, T^ ax s=s Q.QiJ^l 
hence J ~ 156 K. This value is in good agreement 
with our DFT estimate of 150 — 250 K. However, the 
uniform-chain fit overestimates the susceptibility below 
80 K (Fig. 2]). To improve the fit in the low-temperature 
region, several extensions/modifications of the uniform- 
chain model are suggested. 

First, the interchain coupling J' reduces quantum fluc- 
tuations, thus leading to a lower magnetic susceptibility 
at low temperatures. Indeed, we found a good fit of 
the experimental data down to 30 K with J' /J = 0.4, 
J = 150 K, and g = 2.13 (dashed line in Fig. 0J. The 
value of J' is, however, far too large compared to the 
DFT estimate of J' /J ~ 0.01. 

At present, DFT calculations can provide reliable nu- 
merical estimates only for the case of isotropic (Heisen- 
berg) couplings. Therefore, anisotropic effects should be 
considered "on top" of the isotropic, DFT-based model. 
Various simulation techniques (e.g., QMC) are capable 
of including the anisotropy, thereby offering a possibility 
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FIG. 4. (Color online) Magnetic susceptibility (x) of 
(NO)Cu(N03)3 compared to different spin models: the uni- 
form spin chain, coupled uniform chains, the uniform chain 
with Ising anisotropy [Eq. JJJ], and the uniform chain with 
staggered anisotropy [Eq. ((2]), A = 0], see text for details. 
Experimental data are the scaled ESR intensities from Ref. [jj. 



of a direct comparison to experiments. 
The symmetric (Ising/XY) anisotropy^ 



H = J 



<ij> 



( QX QX 



SVSV + (1 + A)S*S*) (1) 



has pronounced effect on the low-temperature part of the 
susceptibility curve. In particular, the Ising anisotropy 
(A > 0) opens a spin gap and reduces the low- 
temperature susceptibility. Thus, we are able to fit the 
data down to 20 K with J = 95 K, A = 1, and g = 2.07 
(solid line in Fig. UJ. Since experimental data for other 
Cu+ 2 compounds suggest A < 0.2 (Refs.[H andSO), the 
fitted value of A = 1 looks rather overestimated. 

Other anisotropy effects are the (/-tensor anisotropy 
and the antisymmetric Dzyaloshinskii-Moriya (DM) ex- 
change. Owing to the P1\jm symmetry of the 
(NO)Cu(N03)3 structure, the neighboring Cu atoms are 
imaged with the 2i screw axis, leading to the possible 
staggered anisotropy of the <?-tensor (according to Ref. 0, 
there is a sizable difference between the g-tensor compo- 
nents: gu = 2.06 and g± — 2.36). The mirror planes are 
perpendicular to the b axis and run between the neighbor- 
ing Cu atoms, thus confining the DM vector to be perpen- 
dicular to the b axis (i.e. to lie within the ac plane). Fur- 
ther on, the 2i screw axis induces opposite DM vectors on 
the neighboring bonds and leads to their staggered con- 
figuration. Following the general framework developed 
in Ref. l4ll we describe the uniform chain with staggered 
anisotropy using the Hamiltonian 



H 



j E s ? 



QX 



QV qV 



(1 + A)S?S* 



-h a ^2sf-h 8 J2(-^s!, (2) 



where the first term is the bilinear exchange with the 
symmetric anisotropy A (i.e., A from Eq. (p} modified 



by the staggered anisotropy). The effective uniform (h u ) 
and staggered (h s ) fields depend on the applied external 
field (H ) and on the staggered anisotropy. 

The general effect of the staggered anisotropy is the 
opening of a spin gap and thus a reduction of the low- 
temperature susceptibility for certain directions of the 
applied field, 41 At fixed A, the gap mainly depends on 
h s . For simplicity, we fix A = and fit the data down 
to 25 K with h s /J = 0.15, J = 160 K, and g = 2.06. 
The uniform component of the field can be varied in a 
wide range, thus leaving freedom for H and the stag- 
gered anisotropy parameters.— Since h s = H sin ^ and 
tana = |D|/J (Ref. EH), we find H/J > 0.13 that defi- 
nitely exceeds the typical field value of 0.3 T in an X-band 
ESR experiment. The |D| value is effective and implicitly 
contains the staggered anisotropy of the g-tensor which 
is presently unknown. 

We conclude that none of the anisotropy parameters 
can be taken as a sole reason for the reduced suscepti- 
bility at low temperatures. However, the combination of 
different anisotropics slightly improves the situation and 
brings the numbers closer to our expectations: A = 0.5 
and h s /J = 0.06 (dashed line in Fig. [5]). Nevertheless, 
the presence of such a strong exchange anisotropy has to 
be challenged experimentally. 

Weak dimerization might be an alternative to the ex- 
change anisotropy. Since exchange couplings are highly 
sensitive to details of the crystal structure, a weak struc- 
tural change will readily open the gap and reduce the 
susceptibility owing to the alternation of the exchange 
couplings along the chain. Experimental data are well 
fitted by an alternating J — J" model^i with J = 176 K, 
the alternation ratio J"/ J = 0.87, and g — 2.10 (solid 
line in Fig. [5]). Although the available structural data 
do not suggest the dimerization, we are not aware of 
any structural studies below 160 K. Moreover, the non- 
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FIG. 5. (Color online) Same as Fig. [4] for the uniform spin 
chain, the anisotropic spin chain [Eq. ((5J], and the alternating 
spin chain. The inset shows bulk susceptibility^ and the fit 
with Eq. © (A=0.5, h s /J=0.06]). See text for details. 



trivial temperature dependence of the ESR linewidth 9 
might also be related to a structural change. 

Trying to get further support for one of the above 
(rather speculative) scenarios, we consider other experi- 
mental data. Specific heat was reported in a narrow tem- 
perature range only (1.8 — 10 K).— It shows a minimum 
around 5 K with an increase towards lower temperatures 
interpreted as a Schottky anomaly (although the char- 
acteristic maximum was not observed) and a typical in- 
crease towards higher temperatures due to the phonon 
contribution which is cubic in T. Since J ~ 150 K, 
the magnetic specific heat will show a maximum around 
70 K)2£ hence the experimental data in this temperature 
range (along with the proper non-magnetic reference to 
estimate the phonon contribution) could be helpful. Re- 
garding the Schottky anomaly, no indication of its intrin- 
sic origin has been given.— 

In contrast to the specific heat, the raw magnetic 
susceptibility data are very instructive. Owing to the 
huge impurity contribution, the susceptibility maximum 
around 100 K is hardly visible. This situation is very 
typical for spin-chain systems, since ID magnets are 
highly sensitive to impurities (a single defect breaks the 
chain). 24 Following previous studies^ we fit the data us- 
ing the expression 



X( T ) = XO + C imp /T + Xchain(r), 



(3) 



where xo accounts for core diamagnetism and van Vleck 
paramagnetism, Ci mp /T is the Curie law to fit the im- 
purity contribution, and Xcha,m(T) is the susceptibil- 
ity of the uniform spin chain* 3 -! We find a remarkably 
good fit down to 2 K with xo = 7.7 x 10 -5 emu/mol, 
Cimp = 0.015 emu K/mol (4 % of spin-i impurities), 
J = 150 K, and g = 2.12. Since the same model poorly 
fits the intrinsic susceptibility from ESR below 80 K (see 
Fig-H]), we are left with two options: i) the impurity con- 
tribution is not properly described with the Curie law 
(then, the origin of this unusual "impurity" contribution 
is worth to unravel) or ii) the ESR bears a systematic 
error that causes the underestimate of x a t low tempera- 
tures (probably, due to the separation of the intrinsic and 
impurity signals in the spectra). To resolve this puzzling 
issue, additional experimental studies, such as suscepti- 
bility measurements on single crystals or nuclear mag- 
netic resonance, are highly desirable. 

Finally, we would like to comment on the possible LRO 
in (NO)Cu(N0 3 ) 3 . Although Ref . H claims the absence of 
the LRO down to 2 K, we suggest a different estimate of 
T^v < 5 K, since the upturn of the specific heat ("Schot- 
tky anomaly" ) may conceal the expected weak transition 
anomaly below 5 K. The weakness of the anomaly is a 
natural consequence of T/v<C./, which results in a signif- 
icantly small amount of entropy released at Xjy Further 
on, ESR intensities diverge below 10 K and might also in- 
dicate the onset of LRO. Our upper estimate of T/v=5 K 
corresponds to T^/J = 0.03. We will show that such 



a low Tjv is typical for a ID system and should not be 
taken as a sole evidence of magnetic frustration. 

In (NO)Cu(N03)3, uniform spin chains with J ~ 
150 K are coupled by J' /J ~ 0.01 and J ± /J ~ 0.0025. 
Unfortunately, the case of spatially anisotropic interchain 
couplings has not been considered theoretically, yet. As- 
suming the same interchain coupling J 1 / J ~ 0.01 along 
the two directions, we arrive at T^ / J — 0.021 (Ref. |44T ) 
which is already below our upper estimate of 0.034^ With 
proper accounting for the spacial anisotropy, T/v should 
be even lower, because it is largely determined by the 
lowest interchain coupling J± as the main obstacle for 
LRO. We conclude that the lack of clear observation of 
LRO in (NO)Cu(N03)3 results from the pronounced one- 
dimensionality of the system. The low T/v does not evi- 
dence the strong frustration. Moreover, the DFT-based 
model approach accounts for all isotropic exchange cou- 
plings and does not show the frustration. 

The low T N in (N0)Cu(N0 3 ) 3 can be compared to 
other spin-i uniform-chain magnets. Weak interchain 
couplings of J'/ J ~ 0.01 were previously observed in 
Sr 2 Cu0 3 (Refs. and \M) and Sr 2 Cu(P0 4 ) 2 (Refs. M 
andl47l). In these compounds, the magnetic ordering tem- 
peratures are Tn/J ~ 2x 10 -3 and 5 x 10~ 4 , respectively, 
thus an ordering temperature of (NO)Cu(N03)3 should 
be quite low and may even lie below 2K (~ 10~ 2 J). 

In summary, DFT calculations suggest a consistent de- 
scription of (NO)Cu(NOs)3 as a uniform-spin-chain sys- 
tem with weak and non-frustrated interchain couplings. 
The magnetic ordering temperature Tn/J is predicted to 
be below 0.03, indicating that LRO could not be observed 
in previous experiments. To find experimental signa- 
tures of the LRO, low-temperature studies and sensitive 
experimental techniques (such as muon spin relaxation) 
should be applied. Regarding the behavior above T/v, 
the bulk magnetic susceptibility data follow the uniform- 
chain model, whereas the intrinsic magnetic susceptibil- 
ity measured by ESR shows lower x below 80 K. We 
note that other spin-chain systems also show puzzling 
behavior at low temperatures. For example, Ref. [48| re- 
ports a complex ESR spectrum of (6MAP)CuCl3 which 
was previously known as a uniform-spin-chain compound. 
To explain the experimental spectrum, the authors of 
Ref. 48 had to consider a spin chain with several inequiv- 
alent exchange couplings, although the structural data 
did not show any signatures of the distortion. It is possi- 
ble that numerous systems, assigned to the uniform-chain 
model from bulk magnetic susceptibility measurements, 
are more complex than they appear. Further studies 
should shed light to this problem. 
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